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The phenomenon of Coherent Population Trapping (CPT) of an atom (or solid state "artificial 
atom"), and the associated effect of Electromagnetically Induced Transparency (EIT), are clear 
demonstrations of quantum interference due to coherence in multi-level quantum systems. We 
report observation of CPT in a superconducting phase qubit by simultaneously driving two coherent 
transitions in a A- type configuration, utilizing the three lowest lying levels of a local minimum of 
a phase qubit. We observe 60 (±7)% suppression of excited state population under conditions of 
CPT resonance. We present data and matching theoretical simulations showing the development 
of CPT in time. Finally, we used the observed time dependence of the excited state population to 
characterize quantum dephasing times of the system. 
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Coherent population trapping (CPT) is a physical ef- 
fect in which a quantum system is driven into a coher- 
ent superposition of two long-lived states by a combina- 
tion of two optical fields. These fields couple each of the 
long-lived states to a common, shorter-lived excited state. 
CPT has been observed for some time in atomic systems 
[l|. CPT refers to the dynamics of the atom while the 
associated effect on the optical fields is referred to as elec- 
tromagnetically induced transparency (EIT) [2] . EIT has 
been shown to give rise to a number of interesting effects, 
including slow light , low light level non-linear optics 
and single photon quantum storage and memory 
41 0| . There have been impressive demonstrations in 
quantum wells 
dots [ill. 
In 



11 . Il2j and more recently in quantum 



14j it was proposed to apply these concepts to the 



microwave, rather than optical, domain in superconduct- 
ing quantum circuits. This step is interesting for two 
reasons. First, it applies CPT and EIT and related non- 
linear optical effects to superconducting qubits, which 
have shown great promise for future quantum informa- 
tion processing systems Second, it also involves 
observing these effects in a completely different region of 
the electromagnetic spectrum. Some experimental steps 
in this direction were taken in [13, [HI , in which the re- 
lated effect of Autler-Townes (AT) splitting [l9| was ob- 
served. 

In this letter, we present direct experimental observa- 
tion of Coherent Population Trapping (CPT) in a three 
level A configuration superconducting system in which 
two long-lived levels (|0) and |1)) are coupled to a short- 
lived higher- lying level (|2)). By resonantly driving two 
transitions of this system with a two-tone microwave 
field, we observe the creation of a coherent superposi- 
tion state of the lower lying levels, which is known as a 
"dark" state and is central to the observation of CPT 
[if. We observe the probe resonance in the absence of 
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FIG. 1: (a) Single well anharmonic potential shown with four 
energy levels. Arrows illustrate the dominant driven transi- 
tions and decoherence channels in our system three level A 
system, (b) Rabi oscillation of the /oi (blue) and /02 (red) 
transitions are shown for potential with /oi = 6.7 GHz. The 
two curves are offset for clarity, (c) A micrograph of the 
measured qubit. The Josephson junctions in the circuit are 
epitaxially grown Rhenium with crystalline AI2O3 tunneling 
barriers, (d) Circuit diagram of qubit, SQUID, and flux bias 
line. The qubit is controlled via the flux bias line and the ca- 
pacitively coupled microwave line. The qubit state is probed 
by current biasing the SQUID and measuring the resulting 
voltage. 



CPT conditions, and a 60(±7)% suppression of this res- 
onant excitation under conditions of CPT. Furthermore, 
we observe the time dynamics of the development of the 
CPT resonance and use the observed time dependence 
to extract decoherence rates, as proposed in [14]. Our 
experimental results are corroborated by theoretical sim- 
ulations of the density matrix evolution of the system. 

In our implementation, the three states of our A sys- 
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tern are the three lowest levels of a metastable phase 
qubit potential well ^] (Fig. [U^a)). This device ex- 
hibits an anharmonic potential in which the transition 
frequencies, /oi and /12, differ by a few percent. The 
design of the qubit seen in Fig. [Hc,d) is identical to 
the devices explored in [2l|. The fabrication of the junc- 
tions was targeted to yield a flux tunable transition fre- 
quency in the 5-10 GHz range. For the CPT exper- 
iment we biased the qubit to have foi = 6.205 GHz 
which makes /12 = 5.865 GHz. This corresponds to a 
shallow qubit potential that supports approximately 8 
metastable bound states. For the sake of visual simplic- 
ity these other levels are not include in Fig. [U^a). We 
choose this flux bias to make the anharmonicity param- 
eter A = (/oi — /i2)/2 = 170 MHz sufficiently large and 
to create a strong damping of |2) which is advantageous 
for observing CPT. 

When a microwave fleld is applied at the frequency /oi 
we obtain Rabi oscillations with a period determined by 
the Rabi frequency Q'^^^^ proportional to the fleld en- 
velope amplitude. To create a A system, we couple |0) 
to 1 2) via a two-photon virtual transition with a probe 
microwave fleld near frequency /02/2 = 6.035 GHz. In 
this virtual excitation scheme the single photon energy 
is signiflcantly detuned from /oi. The detuning is deter- 
mined by the anharmonicity parameter A. The resonant 
Rabi frequency associated with the |1) |2) transition 
is l^p^^^ = V^l^p^^^ and for sufficiently small microwave 
fleld amplitude <C A), the effective two- 

photon coupling amplitude is 

22]. The accuracy of these expressions was conflrmed in 
22] where it was noted that this two-photon virtual tran- 
sition amplitude quantitatively explained the observed 
coupling rates. It was also noted there that the two- 
photon virtual process dominates the direct two-photon 
|0) o 1 2) process. Contrary to the spectroscopic results 
in we observe only the /oi and /02 transitions and 
none of the higher level transitions observed in [22[. In 
our experiment we choose Q^^"^ ^ (27r)48 MHz, giving 
Ap = (27r)6.8 MHz. 

Rabi oscillations shown in Figure 1(b) were measured 
with a qubit bias that exhibited longer coherence times 
[23] than the bias used for the CPT experiment. This 
demonstrates coherence of both the |0) ^ |1) and |0) ^ 
1 2) transitions. The |0) ^ |2) oscillations are slower than 
the |0) ^ |1) oscillations due to the smaller effective Rabi 
frequency, and more damped due to the faster relaxation 
of state 1 2) [23| . We note that in this scheme there is some 
unwanted population excitation of the intermediate level 
\n^p^^\y2A^ ^3%. This process is accounted for in our 
theoretical analysis and it effects our data quantitatively 
but does not prevent CPT from clearly being observed. 

To suppress the probe transition process via CPT, we 
apply a second coupling microwave fleld at frequency /12. 



The amplitude of this transition is given hy Ac = , 

(12) 

where Qc is the resonant Rabi frequency. When both 
flelds are applied simultaneously, the system is ideally 
driven from the ground state |0) into the "dark" state 
superposition lip) = {Ac\0) — Ap\l))/A, where A = 
Y^IApp + l^cP- In this superposition, the excitation am- 
plitudes for each of the individual probe and coupling 
transitions are equal and opposite, resulting in a quan- 
tum interference which suppresses the excitation of state 
|2). 

In order to initialize the system into a well deflned 
state, we let the qubit relax to the state |0). To observe 
CPT we then apply a two-tone microwave fleld at fre- 
quencies fp and /c, with independently controlled power 
at each frequency, and duration to using a pulsing scheme 
similar to [24] . Our microwave pulses were generated by 
sending trapezoidal envelope pulses through a 1 GHz low 
pass fllter and into the IF port of a microwave mixer. 
Analysis of Rabi and Ti data showed that we had inde- 
pendent control over the amplitude and duration of these 
pulses. After the time to, the qubit state is destructively 
measured with a calibrated pulse that probes the pop- 
ulation of the state |2) only (P2) 1251] • We repeat this 
procedure many times at each fp and fc to get statistics, 
and swept these frequencies about the resonances /02/2 
and /12. 

Figure [2ja) presents our results showing the CPT res- 
onance and Fig. [2fc) shows cuts at two values of fc. 
At fc = 5.73 GHz, far from the coupling resonance 
/12, we observe a probe resonance at fp = 6.035 GHz 
with width of approximately 40 MHz. This is primarily 
broadened by the linewidth of the |0) ^ |2) transition, 

702 = l/2Tf + 1/^0°^^. No CPT is observed here as 
the dark state is quickly dephased. A similar feature is 
obtained in the absence of a coupling fleld. By contrast, 
a cut at fc = 5.86 GHz reveals a 60 (±7)% suppression 
of P2 relative to the off-resonant case. The CPT con- 
trast and uncertainty was extracted from a parabolic flt 
of both curves in Fig.[2fc). A salient feature observed in 
Fig. [Sfa) is the diagonal trench of suppressed P2 at posi- 
tions corresponding to the CPT resonance 2fp — fc = foi 
(white dot-dashed line in Fig. [2fa)). 

We compare our experimental results with theoretical 
simulations in which we propagate the density matrix 
equations in time with a Runge-Kutta algorithm. In the 
model we include the flrst four levels of the phase qubit 
potential well and dipole coupling on |0) ^ |1), |1) ^ |2), 
and 1 2) ^ |3) from both applied flelds. The relative 
dipole moments for these transitions are 
The presence of the state |3) was seen to have a negligi- 
ble effect on our simulation results, as was expected. The 
population of state |3) was found to be < 0.5% for all of 
our simulations. The overall calibration between Rabi 
frequency couplings ft^^l^ and microwave powers was flt- 
ted as a free parameter. Since the relative phase of the 
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FIG. 2: Demonstration of CPT in a phase qubit. (a) Popu- 
lation of state 1 2) (P2) measured after a two tone microwave 
pulse with a duration of to = 30 ns and frequencies fp and 
fc. Each population measurement was performed by count- 
ing the number of SQUID switching events out of 4,000 trials 
and accounting for 80% measurement fidelity and 3% back- 
ground counts. The power levels for the probe and coupling 
fields were -6.5 dBm and -15 dBm at the source, respec- 
tively. The white dot-dashed line shows the two CPT res- 
onance 2fp — fc = foi (b) Theoretical data obtained by nu- 
merically solving the four level master equation (see text). 
We obtained the best match with experimental data by set- 
ting the resonant Rabi frequencies for these powers to be 
= (27r)48 MHz, Q^^^^ = (27r)32 MHz. These values 
yield Ap = \nl^^^\^/2A = 6.8 MHz, Ac = ^i^^^ - 32 MHz. 
(c-d) Cross sections of experimental data (with statistical er- 
ror bars) (c) and theoretical data (d) taken at the values of 
fc indicated by the solid and dashed lines in (a) and (b). The 
cuts show the one-photon resonance when fc is off-resonance 
and show the suppression of P2 when fc = 5.86 GHz. 



probe and coupHng field varied randomly from shot to 
shot in our experiment, our theoretical results were aver- 
aged over all possible phase differences of the microwave 
fields. The dipole moment for direct |0) ^ |2) transitions 
is much smaller [22| and furthermore is far detuned from 
both microwave fields. Population decay was included 
according to independent measurements of the energy re- 



laxation times, T{ = 108 ns and T| = 77 ns. We 
extracted the dephasing time T^^^^ = 6(±1) ns from the 
width of several off resonant cuts near fc = 5.73 GHz. 

Results of our theoretical simulations are presented in 
Fig. Efb). We see good agreement with the experimen- 
tal data in Fig.[2fa), with the theory clearly showing the 
broad horizontal stripe corresponding to the probe res- 
onance and the diagonal trough due to CPT. A more 
quantitative comparison is seen by examining the cuts 
in Fig. [2fc) and (d). By fitting data near the CPT res- 
onance at fc = 5.865 GHz, we found best agreement 
for dephasing times T^^^ = 30 ns and T^^^^ = 30 ns. 
Rabi data taken for this potential was consistent with 



this value of T^^\ In the limit of a perfect three- 
level A system and an infinitely long lived coherence 
701 = 1/2T^^^^^^1/T^^^^ = 0, the CPT suppression un- 
der conditions of EIT resonance (2/p — fc = foi) would 
be 100%. The resonant P2 suppression predicted in the 
theoretical simulation is 52%, in fairly good agreement 
with the observation of 60 (±7)%. Analysis of the the- 
oretical results reveals that the finite probability of ex- 
citation is primarily due to dephasing of the |0) ^ |1) 
transition, but also partly due to undesired off-resonant 
coupling [26]. The most dominant such coupling is due 
to the probe field coupling the |1) ^ |2) transition. 

We note that our configuration is distinct from that 
employed in the recent observation of Autler-Townes 
splitting in a superconducting phase qubit [HI, where a 
ladder scheme was used and population was suppressed 
in the intermediate level |1). This leads to two impor- 
tant qualitative differences. First, in contrast to [18|, 
we directly measure the expected population suppression 
(which in our case is of |2)). Second, the fact that we 
create a dark state of the two longest lived states in the 
system allows us to unambiguously observe CPT, rather 
than the closely related phenomena of Autler-Townes 
(AT) splitting pJi,,18]. In AT, the coupling field dresses 
the system, splitting the probe resonance, of width 702, 
into two distinct resonances, spaced by the coupling field 
amplitude Ac ^ 702- In this regime, the suppression 
of P2 at the probe resonance can be attributed to large 
detuning (^ Ac/2) from the dressed resonance peaks. 
By contrast, when Ac < 702 the resonances still largely 
overlap and the suppression of P2 can only be explained 
by a quantum inteference of the excitation processes. 
In our experiment, we estimate Ac = (27r)32 MHz, 
and 702 = 27.5 MHz, which would correspond to only 
1-2702(702+701)/ ((702 +701)' + ^^?) = 17(±13)% sup- 
pression of P2 attributable purely to the AT splitting. 
The uncertainty in this quantity comes from the uncer- 
tainty in the value of T^^\ Using this same model the 
contrast including quantum interference is predicted to 
be 62%, which is consistent with our experimental re- 
sults. 

To further understand CPT and EIT in our system, 
we varied the duration to of the two microwave fields. 
Figure presents measurements of P2 as a function 
of to and fc keeping the probe field on the two-photon 
resonance fp = /02/2. The results show the excitation 
into 1 2) is slowed when the coupling field is at the CPT 
resonance fc = /12. Figure [31^ c-e) shows cuts in fc at 
several different times. Note that the probe is always on, 
resonantly exciting the |0) |2) transition. As time pro- 
gresses the CPT feature narrows and deepens until the 
maximum suppression 55%) is reached at about 20 
ns. After this time, contrast is degraded by the decoher- 
ence rate 701 and in the limit of longer times, the EIT 
notch disappears. 
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FIG. 3: Evolution of EIT as a function of microwave pulse 
duration. In this data set /oi = 6.19 GHz, /12 = 5.85 GHz, 
and /02/2 = 6.02 GHz, 15 MHz shifted from the data set 
in Fig. [21 fp was held fixed at 6.0158 GHz, essentially on 
resonance with the |0) ^ |2) two-photon transition. Each 
population measurement was an average of 5,000 trials, (a) 
Experimental measurement of P2, and (b) theoretical calcu- 
lation of P2. The theoretical resonant Rabi frequencies were 



(01) 



= (27r)48 MHz and Qi^^^ = (27r)35 MHz. (c)-(e) Cross function of T' 



n 

sections of experimental data (symbols) and theory data (line) 
indicated by black dashed lines on (a) and (b). Statistical er- 
ror bars are approximately equal to the symbol size. White 
dot-dashed lines indicate cross sections shown in Fig. [4] 



Since the excitation of P2 is driven by 701, it was pre- 
viously proposed [14] to extract this decoherence rate by 
measuring loss in a CPT experiment. We realized this 
idea by taking a cut of the data in Fig. [3l^a) both off 
{fc = 5.77 GHz) and on {fp = 5.85 GHz) CPT reso- 
nance (Fig. |ll^a,b)). The theoretical data from Fig. [3] are 
then presented as solid black lines. We adjusted the var- 
ious dephasing times and energy level frequencies in the 
theory to best fit the experimental data and found the 



best agreement for 



(01) 



12 ns, T 



(02) 



20 ns , and 



^(12) 



20 ns [27]. 



To understand how sensitive a probe this measurement 
is of T^^^ we performed several simulations in which we 
varied this parameter. The inset in Fig. IH^b) shows P2 
at to = 40 ns as a function of Tj°^^ We see that in 
the regime (4 ns < T^^^ < 30 ns) we are able to mea- 
sure T^^\ At longer T^^^ the population asymptotically 
approached a level P2 ~ 0.05 which is dominated by off- 
resonant coupling effects [26j . Performing our experiment 
in a regime with greater anharmonicity would suppress 
this effect, allowing measurement of longer T^^^ and, 
in the event of obtaining larger coherence times, greater 
CPT suppression of P2. 

In summary, we have directly observed CPT in a three- 
level superconducting phase qubit and seen suppression 
of the excited level population that is clearly attributable 



FIG. 4: Cuts in time of the experimental data (with statistical 
error bars) presented in Fig. [^indicated by white dot-dashed 
lines in a,b) at two different frequencies fc corresponding to 
(a) off CPT resonance /12 - fc = SO MHz and (b) on CPT 
resonance f 12 — fc = MHz. Also shown are theoretical 
simulations of the evolution of P2 over time shown for three 
different values of T^^'^^ with other parameters the same as 
in Fig. [3l (inset) P2 for fc = /12 for to = 40 ns plotted as a 
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The black dashed line shows the population 
= 80 MHz). 



to quantum interference. We did this by creating a A 
configuration using a two-photon virtual probe transi- 
tion and monitoring the population of |2). We have 
matched our results with theoretical simulations of the 
full density matrix equations, including expected deco- 
herence and cross-talk effects. We have furthermore ob- 
served the development of the characteristic EIT notch 
in time and used these time dynamics to characterize the 
decoherence time of the |0) ^ |1) transition. 

These results are of fundamental interest as a clear 
observation of a quantum interference at a macroscopic 
scale. Furthermore, they are directly pertinent to future 
investigations of low light level nonlinear optics in the mi- 
crowave domain and quantum information applications 
utilizing superconducting quantum circuits. 
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